Antiferromagnetism of the 2D Hubbard Model at Half Filling: Analytic 

Ground State at Weak Coupling 

f^ ' Michele Cini and Gianluca Stefanucci 

^^ ' Istituto Nazionale di Fisica della Materia, Dipartimento di Fisica, 

<«_^ , Universita' di Roma Tor Vergata, Via della Ricerca Scientifica, 1-00133 
Cn ' Roma, Italy 



Ph. 






c^ 



o 



^ 

c 



X 



c^ 



We introduce a local formalism to deal with the Hubbard model on a, N x N square 

lattice (for even A^) in terms of eigenstates of number operators, having well defined 

point symmetry. For U -^ 0, the low lying shells of the kinetic energy are filled in the 

"^J I ground state. At half filling, using the 2N — 2 one-body states of the partially occupied 

/ 2A'' — 2 \ ^ 
shell Shf, we build a set of I , , I degenerate unperturbed ground states with 

Sz = which are then resolved by the Hubbard interaction W — U'^ hrihri- In Shf 

we study the many-body eigenstates of the kinetic energy with vanishing eigenvalue 

. ■ of the Hubbard repulsion {W = states). In the Sz = sector, this is a A'' times 

degenerate multiplet. From the singlet component one obtains the ground state of the 

Hubbard model for U — 0'^ , which is unique in agreement with a theorem by Lieb. 

The wave function demonstrates an antiferromagnetic order, a lattice step translation 

T-r^ ' being equivalent to a spin flip. We show that the total momentum vanishes, while the 

pH , point symmetry is s or d for even or odd N/2, respectively. 

o 
o 

I. INTRODUCTION 

^ ■ The discovery of high temperature superconductors ||l| enhanced the interest in models of two- 

OO I dimensional, strongly correlated electron systems, such as the 2D Hubbard model. The (repulsive) 

Hubbard hamiltonian is highly idealised, for istance two particles are allowed to interact only on-site, 

yet the model already displays an interesting phase diagram. The fluctuation exchange (FLEX) g 

^-*) ' . diagrammatic approach, which is based on a conserving approximation and, independently, renormal- 

(^ ' ization group techniques H M show that near the antiferromagnetic phase at half flUing, there exists 

f^ , a superconducting phase too, with a momentum dependent gap, of d wave symmetry. Indication 

of a possible instability of the Fermi liquid towards pairing near half filling also comes from cluster 

2 ! diagonalizations 0- Therefore, exact results on the half filled Hubbard model may be relevant 

to antiferromagnetism and to the mechanism of the superconducting instability as well. 

In the strong coupling limit the double occupation of the same site is energetically suppressed and 
the model at half filling is equivalent to the Heisenberg model with an antiferromagnetic exchange 
Q ' interaction |g]. A popular approach takes care of the strong repulsion between two opposite spin 

Q , fermions by a Gutzwiller Q projection, i.e. by throwing out of the Hilbert space the double occupation 

states. 

However, truncating the Hilbert space in this way costs lots of kinetic energy, so at finite U the 
system must allow double occupation, also in the ground state. At weak coupling it makes sense 
^ I to speak about particles in filled shells, which behave much as core electrons in atomic physics, and 

particles in partially filled, or valence, shells. Remarkably, particles in partially filled shells can totally 
avoid double occupation at no cost in energy; they do so, forming W=0 states, that are defined as 
many-particle eigenstates of the kinetic energy with no double occupation. Below, using a 
new formalism, we show how W = states arise by symmetry. 

An important theorem on the Hubbard model at half filling is due to Lieb |1Q] : the ground state 
for a bipartite lattice is unique and has spin ^||5i| — |52|| where |iSi| {\S2\) is the number of sites in 
the Si {S2) sublattice; here and in the following, |5| will be the number of elements in the set S. It 
is worth to observe that the theorem makes no assumptions about the symmetry of the lattice. For a 
N X N square lattice, with N even, the ground state is a singlet and in Ref. flit] it was shown that in 
the strong coupling limit it has total momentum Ktot = (0, 0) and s wave (a;^ +2/^) or d wave (x^ — y'^) 
symmetry for even or odd N/2 respectively. 

In this paper we build the exact ground state of the Hubbard model at half filling and weak coupling. 
In Section |B we state more precisely the problem we want to solve and we define some notations. In 



Section II] we show how for each site we can build a local one-body basis set, which is well suited to 
write the antiferromagnetic many-body wave function. Finally in Section |^ we explore the symmetry 
properties of the ground state, and deduce the same quantum numbers as predicted by Refs. pO| | pTl 
at strong coupling. 

The present approach lends itself to obtain exact results for other fillings as well, but this will be 
shown elsewhere 



II. THE HUBBARD MODEL AT HALF FILLING AND WEAK COUPLING 

Let us consider the Hubbard model with hamiltonian 

a {r.r') r 

on a square lattice of N x N sites with periodic boundary conditions and even N. Here a =t, i is 
the spin and r, r' the spatial degrees of freedom of the creation and annihilation operators c^ and 
c respectively. The sum on (r, r') is over the pairs of nearest neighbors sites and fira is the number 
operator on the site r of spin a. The point symmetry is C4y, the Group of a square ||I^; besides, 
H is invariant under the commutative Group of Translations T and hence the Space Group |Q 
G = T (g) Civ', ^ means the semidirect product. We represent sites by r = {ix, iy) and wave vectors 
by fc = (kx, ky) = ^{ix, iy), with ix,iy = 0, . . . , N — 1. In terms of the Fourier expanded fermion 
operators Cka = jj- X^r ^^'^^^rcr, we have Hq = ^j, e(fc)c|,g.CfcCT with e{k) — 2t[cos kx + cos ky]. Then the 
one-body plane wave state c|,^|0) = {ka) is an eigenstate of Hq. 

We study the ground state. Let Shf denote the set (or shell) of the k wave vectors such that 
e{k) — 0. At half filling {N'^ particles) for U = the Shf shell is half occupied, while all |fc) orbitals 
such that e(fc) < are filled. The k vectors of Shf lie on the square having vertices (±7r, 0) and 
(0, ±7r); one readily realizes that the dimension of the set Shf, is \Shf\ = 2N — 2. Since N is even and 
H commutes with the total spin operators, 



S, 



i^(fM-n,i), S+=J24-^Cri, S-^iS+)\ (2) 



at half filling every ground state of Hq is represented in the Sz — subspace. Thus, Hq has 

IN -2\'^ 

\j , ) degenerate unperturbed ground state configurations with Sz = 0. We wish to study 

below how this degeneracy is removed by the Coulomb interaction W already in first-order pertur- 
bation theory. Actually most of the degeneracy is removed in first-order, and with the help of Lieb's 
theorem we shall be able to single out the true, unique ground state of H. In Appendix |A| we show 
that the structure of the first-order wave functions is gained by diagonalizing W in the truncated 
Hilbert space H spanned by the states of A^ — 1 holes of each spin in Shf- In other terms, 
one solves a 2N — 2-particle problem in the truncated Hilbert space Ti. and then, understanding the 
particles in the filled shells, obtains the first-order eigenfunctions of H in the full A^^-particle problem. 
We underline that the matrix of Hq in Ti. is null, since by construction Ti. is contained in the kernel of 

Hq. 

The operator J^r nr^nri has eigenvalues 0,1,2,... and so the lowest eigenvalue of W is zero (in 
other terms, W is positive semi-definite). The unique ground state of the Hubbard Hamiltonian for 
[/ = 0+ at half filling will turn out to be a M^ = singlet state of 2A^ — 2 holes in Shf (filled shells 
being understood). We shall obtain the W = Q states G H. It is clear that, although the C/ = case 
is trivial, at C/ = 0"*" we are still facing a bona fide many-body problem, that we are solving exactly 

ra. 



HI. LOCAL FORMALISM FOR THE GROUND STATE WAVE FUNCTION 

In the present section we first define a basis of local orbitals; then, we demonstrate a method 
for actually constructing the basis in the general case. We use the 4x4 case as a simple example 



and then generalize. Then we show that using the local basis, the many body wave function of the 
antiferromagnetic ground state can be projected out as the singlet component of a single determinant, 
which is amazingly simple for an interacting system. 



A. The Basis Set: Definition 

Since W depends on the occupation number operators fir , it is intuitive that its properties in Tl are 
best discussed by a suitable one-body basis of Shf such that at least one of these operators is diagonal. 
In addition, a convenient basis should exploit the large G symmetry of the system. If Shf were a 
complete set (A''^ states), one would trivially go from plane waves to atomic orbitals by a Fourier 
transformation; instead, we must define the local counterparts of plane-wave states using only the 
2N — 2 states that belong to Shf- For each site r we diagonalize the number operator n^; moreover, 
since fir is compatible with the operations of the point symmetry group C^v we also diagonalize the 
Dirac characters of the Group. The set of Dirac characters defines the irreducible representation 
(irrep); thus we write the one-body basis states {lipQ )} where a comprises the fir eigenvalue and a 
C4i, irrep label. However, different sites yield different sets; the eigenvectors |</?q ) of 71^=0 and those 
\ipa ) of other sites r are connected by unitary transformations. Introducing the primitive translations 
of the lattice ex — (1, 0) (one step towards the right) and Cy = (0, 1) (one step upwards) the primitive 
unitary transformations read: 

2N~2 2N-2 

(3=1 13=1 

The translation matrix T/ knows all the G symmetry of the system, and will turn out to be very 
special. Using such a basis set for the half filled shell the antiferromagnetic order of the ground state 
comes out in a clear and transparent manner. 



B. Teclinique for Building the Basis Set 

To accomplish that, it is not actually necessary to diagonalize any (2iV — 2) x (2iV — 2) matrices. 
The number operator fir = c].Cr (for the moment we omit the spin index) is dealt with most easily by 
the following 

Theorem: Let S be an arbitrary set of plane-wave eigenstates {\ki)} of Hq and (jir)ij — 
{ki\hr\kj) = -|y2 e*'*^' '^^''^ the matrix of fir in S. This matrix has eigenvalues Ai — jjj and 
X2 = ■ ■ ■ — A|5| = . 

Note that \S\ < N"^; if \S\ = N"^ the set is complete, like the set of all orbitals, and the theorem is 
trivial (a particle sitting on site r is the rir eigenvector with eigenvalue 1); otherwise the theorem is 
an immediate consequence of the fact that (see Appendix H) 

dei|(n.),,-A^,,| = (-A)l'5|"i(l^-A), Vr. (4) 

It is easy to verify that for r = the eigenvector with nonzero eigenvalue is just the totally symmetric 
superposition of all the {|fci)} G S. 

Next, the large set Shf breaks into small pieces if we take full advantage of the G symmetry. Any 
plane-wave state k belongs to a one-dimensional irrep of T; moreover, it also belongs to a star of k 
vectors connected by operations of 6*4^, and one member of the star has k^ > ky > 0. We recall that 
any k e Shf lies on a square with vertices on the axes at the Brillouin zone boundaries. Choosing an 
arbitrary k e Shf with k^ > ky > 0, hence k^ + ky = tt, the set of vectors Rik G Shf-, where Ri € C^y, 
is a basis for an irrep of G. The high symmetry vectors fcyi = (7r,0) and ks = (0,7r) are the basis of 
the only two-dimensional irrep of G, which exists for any N. If A^/2 is even, one also finds the high 
symmetry wavevectors k = (±7r/2,±7r/2) which mix among themselves and yield a four-dimensional 
irrep. In general, when k is not in a special symmetry direction, the vectors Rik are all different, so 
all the other irreps of G have dimension 8, the number of operations of the point Group 64^. 



Below, we shall need the number of these irreps. Since 8 times the number of eight-dimensional 
irreps + 4 times that of four-dimensional ones + 2 for the only two-dimensional irrep must yield 
\Shf\ — 2iV — 2, one finds that S^f contains iVg = i(-^ — 2) irreps of dimension 8 if N/2 is even and 
No — ^{y ~ ^) irreps of dimension 8 if N/2 is odd. 

We note incidentally that G cannot explain the degeneracy 2 A^ — 2 of Shf , because the maximum 
dimension of its irreps is 8. Indeed, the accidental degeneracy of several irreps is due to the presence 
of extra symmetry, i.e. G is a subgroup of the Optimal Group defined in Ref. Ilfl]. 

In this way, Shf is seen to be the union of disjoint bases of different irreps of the Space Group. This 
break-up of Shf enables us to carry on the analysis and build the basis for any N. For illustration, 
we will first consider the case N — A and then generalize. 



C. Example: Basis Set and Ground State for the 4x4 Square Lattice 



As already noted, fc^ = (tt, 0) and fc^ = (0, tt) belong to Shf and are the basis of a two-dimensional 
irrep of G. The 2x2 matrix {nr=o)ij — {ki\nr=o\kj) , with i,j — A,B, has the eigenvector 
\iIj^ ') — -y=(\kA) + Iks)) with eigenvalue Ai — 1/8, in agreement with the above theorem. The 

second eigenvector, with vanishing eigenvalue, is \ipg ) — -y=[\kA) ~ I^b))- As the notation implies, 
both are simultaneously eigenvectors of the Dirac characters and carry symmetry labels; actually the 
symmetries Ai and Bi could have been predicted without diagonalization because the two-dimensional 
irrep of G breaks into Ai © Bi in C^v ■ 

Translating by r, plane wave states pick up a phase factor: |fc) -^ e*'^''|fc). Thus, the eigenstates of 



Ur are \ip^\ 



^{e^k^-r\kA) 



^iksT 



l^s)) with Ai 



with A2 = 0. The primitive translations (H) are performed by 



1/8 and It/;^^'-') = J^{e^^^'r\kA) 



V2^ 



^IfcB)) 



\^i 



J=Ai,Bi 



E 

J=Ai,B_ 



iv/;*°^>(T0,7/, 



x,y 



(5) 



with I = Al, Bi. Using eq. (|5|) one finds the antidiagonal translation matrices T; 



Ty = 



1 

1 



(6) 



So, the orbital of Ai symmetry at r = has Bi symmetry around the nearest neighbour sites, and 
conversely. In particular, \ipj^ ) has vanishing amplitude on a sublattice and \i/jg ) on the other. 



The two-body state \ipj^ /o-IV'b )-o- has occupation for spin a but not for spin — cr on the site r 
under a lattice step translation it flips the spin and picks up a (-1) phase factor: 



0; 



/,"(0)\ 



'(0) 



Ha\')'^\K')-^ ^ \^b:')'^\K')-^ = -ivv:o-.i^sro.; 



(7) 



therefore it has double occupation nowhere and is a VF = state (more precisely, a. W — pair M^ 

For TV = 4, Shf also comprises the basis fci = (7r/2, 7r/2), k2 = (— 7r/2, 7r/2), ^3 = (7r/2, — 7r/2), fc4 = 
(— 7r/2, — 7r/2) of the 4-dimensional irrep of G. This irrep breaks into Ai ® B2 ® -E- in C4V, and such are 
the symmetry labels of the eigenvectors of nr-=o- We easily obtain them using the projection operators 
of C4V. Letting / = 1,2,3,4 for the irreps Ai, B2, E^, Ey respectively, we can write down all the 

eigenvectors of {ki\nr=o\kj) ., with i,j = 1, ... ,4, as j-;/'/ ) = Yl,i=i 0'ii\ki), where O' is the following 
4x4 orthogonal matrix 



O' 



1111 

1-1-1 1 

1-1 1-1 

-1-1 1 1 



(8) 



The state with non- vanishing eigenvalue is again of Ai symmetry. Translating by r we readily get the 

' (r) 

eigenstates |V'/ 

m)ji 



i^T'i^r' 



of rir and of the Dirac characters; in this way, we calculate the translation matrices 
which are block-antidiagonal: 



T^ = 



i 

~i 

i 

-i 



T„- 



-i 

i 

i 

-i 



(9) 



These 4x4 translation matrices are again very special; they are such that for each lattice step 
the subspace of Ai and B2 symmetry is exchanged with the one of E,j: and Ey symmetry, and this 

means that we are about to obtain new W = states. Indeed, jV'yJ ip^ /o'I^b ^e )~cr is a 4-body 

eigenstate of W with vanishing eigenvalue: under a lattice step translation this state does not change 
its spatial distribution but a — > —a. We may write that 



I^A,V'k)-l^k^i^ 



WEjEj^^A.^k)--^ = Wm^bJ-^I^'eJ'eJ^ 



(10) 



and since ip'j^ is the only orbital having amplitude at r = it is evident that this state has double 
occupation nowhere. 

At this stage we have distinct local bases for both irreps of G having their basis vectors in Shf ■ From 

eqs.(||)(Q) we see that under a lattice step translation the subspace spanned by {|V'a /'I^a )' \^b )} 

is mapped in the one spanned by {\4'b )j I'^'e )' \^e )} ^^'^ conversely. Now we are in position to get 
the exact ground state of the 4x4 square lattice at weak coupling. Consider the 6-body determinantal 
eigenstate of the kinetic term Hq 



I'^A} 



\i^7yi:'^^'i:\H7yZ'i^T) 



(11) 



In this state there is partial occupation of site r = with spin cr, but no double occupation. A shift 
by a lattice step produces the transformation 



\<^A} 



l^AI 



(12) 



that is, a lattice step is equivalent to a spin flip, a feature that we call antiferromagnetic property. Since 
the spin-flipped state is also free of double occupation, \^af)<t is a VF = eigenstate, and belongs 
to the first-order ground-state multiplet. Moreover, the single deteminant with the antiferromagnetic 
property may be analysed in its spin components, which must likewise be free of double occupation. 
We show below that there is at least one W = state in H for each S. By Lieb's theorem, the unique 
ground state of the Hubbard model is the singlet component of (|ll|); we shall deal with the projection 
in Sect.fvl. 

0, 



We note that jf/'^ ) and |?/;^ ) are two one-particle states having nonvanishing occupation at r 



and we are getting a new one for each irrep of G; for some applications we might prefer having only one 
such state in Shf- This is easily accomplished after this preparation; shall call \ip) the new local basis. 
According to the above theorem, fir has eigenvalues 3/8 and (5 times) 0. For r = the eigenvector of 
occupation 3/8 is just the totally symmetric superposition of all the |fc) states in Shf, in terms of the 
eigenstates introduced above for the two-dimensional and four-dimensional irreps of G, we may write 



it in the form \ip\ 



(0)x 



"(0)\ 



TSI^A. 



\t 



(0)\ 



The above theorem also grants that the orthogonal 



.(") 



(0)\ 



We 



linear combination of Ai eigenstates of n^^Oi \'P2 ) = y il'^A ) ^IV'a ) ha-s eigenvalue. 

are finished with the ip basis, because the remaining local states are just eigenvectors with vanishing 
eigenvalue, and we may set |(/9^°^) = I'lp^^''), \ipf^) = iV's^"'), lvi°^) = lV'i°^) and |(/7^°^) = IV^i"^). As 



(0)\ 



anticipated, \ip) local basis has the advantage that \ipl ) is the only element with nonzero occupation 
at r = 0; it also preserves the other useful properties that we have analyzed above. So the 6x6 



translation matrix (T; 



l)ij 



n 






1 

VI 

2 
3 




'V3 





i^p 


4' 


) is a 


- 


1 


^^|i 


- 



p2 
V 3 




"73 






















-i 









is a block antidiagonal matrix: 





-i 




T„ 









73O 








\/i 








i 


1 

73 


\/l 








i 





3 


73 






V3 








(13) 



„W\ 



J^ 



{\^',>)A^'2')A^r)h and 5, 



Mx 



The |(p) local basis at any site r splits into the subsets Sa 
{|(/?4 '), \ipi[ ), Iv^g )}; a shift by a lattice step sends members of Sa into linear combinations of the 
members of 5b, and conversely. For the present 4x4 case, we could have obtained the (p basis some- 
what more simply by direct diagonalization on the whole set Shf, but the present approach has the 
advantage of being viable at large N. 

Indeed, \ipl V2 ) i^ equivalent to Itjj^ ^^ ), because this is just a unitary transformation of the 
Ai wave functions. Thus, we may write: 



1^ \ I (0) (0) (0)> I (0) (0) (0), 

\'^AF)a = In V2 V3 )<rm V5 Ve ) 



(14) 



Besides being useful for the sake of illustration because of its relative simplicity, the 4x4 case 
can be thorougly explored on the computer, since the size of TC at half filling is 400. We have used 
Mathematica to diagonalize H + ^S^ , where a small ^ is a numerical device to keep the different spin 
components of the ground state separated. In this way, we observed the fourfold degenerate, W = 
ground state which ^ separates into its singlet, triplet, quinted and septet components, as expected, 
with the separation growing like U^. The antiferromagnetic property of the wave functions was also 
easily and nicely borne out by the numerical results. 



D. The N X N Square Lattice for general Even N 



As discussed in Sect. [II B , we break Shf in the bases of irreps of G it contains. Each basis consists of 
plane- wave eigenstates {|fci)} of iJo and is converted in a Zoca/ one-body basis at site r by diagonalizing 
fir and Dirac's characters . For A^ > 4, Shf contains k vectors that do not possess any special symmetry 
and we get eight-dimensional irreps of G since Rik are all different for all Ri S C^y. In other terms, 
any eight-dimensional irrcp of G is the regular representation of C^y. Thus, by the Burnside theorem, 
it breaks into Ai (B A2 (B Bi (B B2 ® E (B E , with the two-dimensional irrep occurring twice; these are 
the symmetry labels of the local orbitals we are looking for. Now let k G Shf- The only eigenvector 
of the matrix {Rik\hr=o\Rjk) corresponding to the non- vanishing eigenvalue A belongs to Ai. Let 
Ri, i = 1, . . . ,8 denote respectively the identity 1, the counterclockwise and clockwise 90 degrees 
rotation C4 , C4 , the 180 degrees rotation C2, the reflection with respect to the y = and x = 
axis ax, <Jy and the reflection with respect to the x = y and x = —y diagonals a'j^, &_. We can write 

down the eigenvectors of the above nr=o matrix as \ip\ ) = 'Ylii^i Oii\Rik), where k^ > ky > and O 
is the 8x8 orthogonal matrix 



O 



1 

71 



-1 


-1 




-1 


-1 








-1 


-1 




-1 


-1 




— 1 




— 1 


— 1 




— 1 










— 1 


— 1 


— 1 


— 1 


— 1 


— 1 








— 1 


— 1 


-1 




-1 




-1 




-1 


— 1 








— 1 


— 1 





(15) 



is the Ai, B2, Ex, Ey, A2, 



Here, denoting by E' the second occourrence of the irrep E, I — 1, . . 

Bi, E',j,, Ey irrep respectively. A translation by r yields the eigenstates \ipj ) of {Rik\hr\Rjk) and of 



the Dirac characters of the point symmetry group. After very long but elementary algebra one finds 



that the translation matrices {Ti)jj — {ipj \ipj ) are 



T, = 





cosk^ 

0.. 

isink^ 

cos kx i sin k^ 

cos kx i sin k^ 

i sin kx cos kx 

— z sin kx — cos k^ 



coskx 


i sin kx 



i sin kx 

—i sin kx 



fXirf COfe I^X 

































- cos kx 



(16) 



and 



Ty^ 













— cos kx 



— i sinfcr 





I sin kx 





~ coskx —ismkx 

— cos kx Q i sin kx 

i sin kx Q — cos kx 

i sin fca; cos kx 



isinfc^ 

-cosfc:^ 

cosfc^r 








cos kx i sin kx 







(17) 



where we have taken into account that at half filling ky — n — kx- As for the two and four-dimensional 
irreps the translation matrices are in an antidiagonal block form; in particular they are such that the 
set containing the irreps Ai, B2, E is mapped in the set containing the irreps A2, Bi, E' and vice 
versa. This means that if we put 4 particles with spin a in the former 4 irreps and 4 particles with 
spin — cr in the latter 4 ones we obtain an 8-body state for which the translation by a lattice step is 
exactly equivalent to a spin-flip (antiferromagnetic property): 



l<Vi°>i°i<Vl<'<>g<^) 



l<i<>l^A^i^).l<V^"i<V^ 



(18) 



Now recall that the occupation number vanishes for all the local states except the one of symmetry 
Ai. Since in each site the Ai state of spin a does not have the partner of the same symmetry with 
spin —a, the 8-body state of eq.(nq) cannot have double occupancy on any site and therefore it is a 
W = state. 

From now on we shall be engaged with the explicit construction of many-body W — states at half 
filling and it will be understood that we are using the local basis of the site r = 0, so that \ip'-^'>) = \ip). 
We have observed above how the antiferromagnetic property of a many-particle determinant ensures 
that it is a M^ = state. In this way, we easily obtain a determinantal ground state of W in Ti,, i.e., 
at half filling, by creating holes in all the local orbitals of all the irreps, half with spin up and half 

with spin down. Let I^A/ ) be the one-body eigenstate of nr=o belonging to the irrep / of C^y, in 
the space spanned by the basis functions of the rn-th eight-dimensional irrep of G. For even A^/2,the 
four-dimensional representation of G exists and the W — state wave function for the half filled case 



N^ 



N, 



1$^^). ^ i(n <cvL:vL7)^:.>k^i)^i(n ^i.<vgvr)'^kV'k.^Bj-., 



(19) 



with cr =t, i belongs to the first-order ground state multiplet (filled shells are understood, of course) 
For odd N/2, on the other hand. 



No 



Na 



\^af). ^ i(n <<vtvL7)V':i,>.i(n ^tvL":vgvL7)vA)-- 



(20) 



m— 1 



We can see from eqs.(|^)(|To|)(E8|) that |$af)o- flips the spin and picks up a phase factor (see below) 
for each lattice step translation. Therefore, it manifestly shows an antiferromagnetic order [antifer- 
romagnetic property). These results generalize Equation (lO). In both cases, the —a orbitals belong 



to nr=o = 0, and the antiferromagnetic property grants that the states are W = 0. Equivalently, we 
could have obtained a generahsed version of (|lj) by building a basis of symmetry adapted eigenvectors 
of hr on the whole set S^f, which can be done without handling large matrices. 

A few further remarks about |$af)(t a-re in order. 1) Introducing the projection operator Ps 
on the spin 5* subspace, one finds that Ps\^AF)a = \^af)<^ ¥" 0, VS* — 0, ... ,iV — 1. Then, 
a{^AF\W\^AF)a = J2s=i <y (^afI^I^af) <y = 0, and this implies that there is at least one W ^ 
state oi W in H for each S. The ground state of H at weak coupling is the singlet \^'\p)a- 2) The 
existence of this singlet W — ground state is also a direct consequence of the Lieb theorem ||l^ . 
Indeed the maximum spin state \^^p^)„ is trivially in the kernel of W; since the ground state must 
be a singlet it should be an eigenvector of W with vanishing eigenvalue. 3) The above results and 
Lieb's theorem imply that second and higher order effects split the ground state multiplet of H and 
the singlet is lowest. 4) The Lieb theorem makes no assumptions concerning the lattice structure; 
adding the ingredient of the G symmetry we are able to explicitly display the wave function at weak 
coupling. 

In the next section we study the symmetries properties of the singlet component of \^AF)a- 



IV. SPIN PROJECTION AND SYMMETRIES OF $4f> 

The W — state \^AF)cr is a 2(A^— l)-body determinantal state with Sz — but is not an eigenstate 
of the total spin operator S*^ . The various spin components are degenerate in first-order perturbation 
theory, but when higher-order effects are allowed the singlet component is lowest; if we wish to study 
the ground state of the Hubbard model we must project on the singlet. The spin projection operators 
Ps are well known and are reviewed in Appendix O for the sake of clarity. In order to find out the 
good quantum numbers of the ground state, the antisymmetric form 

|. , _ \<^AF)a-\<^AF)-a , . 

I'^AF) = ^= (21) 

is more convenient to work with than the single determinant \^AF)a] using the explicit form of Ps=o 
one finds that the projection is the same. 



A. symmetry under translations 

Eqs. ([7|)([lO| )([lq) tell us that under a lattice step translation each of the determinantal states |$a_f)ct 
of eqs.( p9|j20p undergoes a spin flip, which does not change the irreps of C4V but modifies the order in 
which they appear in the many-body state. Since the fermion operators anticommute, the translated 
determinant is (— 1)^^^|<I>a_f)-(t = — |<&a_f)-o-; but in view of Equation (pi|), |$a_f) -^ \^af) under 
a lattice step translation. Thus |<I>a_f) is an eigenstate of the total momentum with eigenvalue Ktot = 
(0,0). Since the spin projection cannot change this quantum number, it holds for I'l'^'^^) too. 



B. reflections and rotations 

Now we study how |$a_f) transforms under reflections and rotations with respect to the center of 
an arbitrary plaquette of the square lattice. We are not compelled to refer the operations to the center 
of a plaquette, rather than to a site, to characterize the symmetry properties of j^^ii?); indeed the 
system is C^y invariant in both cases. The only reason is to make contact with Ref. pJ^ . 

Since we represent sites by r = {ix,iy) with ix,iy — 0, . . . , A^ — 1, we may choose the center at 

Tpiaq ~ (1/2, 1/2). Let R^"''' denote the C4„ operations with respect to rpiaq and Rt the ones with 
respect to the origin (0, 0). Then for every vector r of our lattice we have from elementary geometry 

Rf'"'r^R,{r-rpiaq) + rpiaq- (22) 



This implies the transformation law for plane- wave states |fc): 

\k) = ^ 5] e-^'^H'') ^ K'"'fc) = ^^e-^'^n^f'^'O = e-''=('--'--«»^''-'"-)|i?,A:), (23) 

r r 

where the last equality can be obtained with a change of variables. By means of eq.(E3|) it is possible 
to know how each irrep of the space group G transforms. 

Two-dimensional irrep: Let us first consider the two-dimensional irrep whose basis vectors are 

|V'Ai> = 75(1^^) + l^s>) and |^^^> = ^{\kA) - Iks)). Under Rf""" we have 

IV^J - |i?r"V^,) = ^(K'"^fcA) - iRf-^ks)). (24) 

The transformed local states |i?^ "''V'/), with / — Ai,Bi, can be expressed in terms of the original 
ones \ip'l): 

i<"^v;')= E mwiRf^'^i) (25) 

J=Ai,Bi 

Using eqs.(p3|) and computing the overlaps {ipj\R^ "''V'/) of eq.(25) we have studied how the W — 
pair state in eq.(R) transforms under a 90 degrees rotation C4 and a reflection with respect to the 
y — and x = y axis ct^; and a'_^_ respectively. After some algebra it is possible to show that 



(26) 



where < — > means that the left hand side transforms in the right hand side and conversely under the 
operation specified on the right. 

For a 2 X 2 square lattice the two-dimensional irrep is the only one in Shf and |$a_f) is explicitly 
given by 

\^af) = mj.l^'I^J-. + WI^J'^Wm)-^- (27) 

By means of cqs.(|2q) is not hard to see that I^&af) transforms as a d wave of x"^ — y^ symmetry. This 
symmetry property cannot change after the spin projection of I^^if) on the S = subspace and we 
conclude that |<&^]^°) transforms as a d wave too. Finally we observe the the 2x2 case is special 
because |$af) is already a singlet and then must coincide with |$^^"). This will be no more true if 
N>4. 

Four-dimensional irrep: Similarly we can study the behaviour of the four-body W = state of 



WJ.IV^J-.- 


- -Wk)'r\^M)-^ 


O4 


M,>.l^^,)-.- 




CTx 


WV^\nj-<r^ 


^-WkJMk)-^ 


<■ 



eq.(lO) under the three operations in eq.(|26|). One finds: 

I^A,^k)-l^k^k)- "^ \^'eJ'eJ<^\^m^'bJ-^ <^^ (28) 

A remarkable feature follows from the above transformation properties. In the A^ = 4 case Shf 
contains the irrep of dimension 2 and the one of dimension 4 so that |$^i?) is 

\^Af) = IV^AiV'B^V'Aj^lV'kV'F^V'Bi)-^ + IV'k^F„V'Bi)<T|V'^iV'B2V'Ai)-<T- (29) 



Using eqs.(p6|)(|2^) it can be shown that |<i>AF) transforms as an s wave of x^ + y^ symmetry for 
iV = 4. Therefore the symmetry of |$a_f) depends on N and the next step will be to determine how 
this happens. To this end we need the transformation properties of an arbitrary 8-dimensional irrcp 
of G contained in Shf- 

Eight-dimensional irrep: 

After very long but simple algebra we have found that the eight-body state of eq. (|8|) transforms as 

|V'Ai'0S2V'Bx'0Bj<T|V'A2'0SiV'£;iV'£;;)-<T < > \lpA2'>pBi'>pE'^'>pE0a\i'Aii^B2->JjE^i'Ey) -a C^'^^ 

\tpAitpB2tpE^tpEy)a\lpA2i'Bi1pE'yE')~a < > \lpA2i^Bill^E'^1pE')a\tpAi1pB2lpE^1pEy) -a CTx 

(30) 

\lpAi1pB2-lpE^tpEy)a\i^A2'ipBi1pE'yE'J~a < > |-0Ai V'Bs V'£;x'«/'isJ<T IV'^a V'SiV'S^ V'B;)-^ 0"+ • 

For kx — > 7r/2 the 8-dimensional irrep is equivalent to two 4-dimensional ones and the eqs.(pO|) are the 
"square" of the eqs.(Eq). Analogously for fc^; ^- tt we obtain four times the 2-dimensional irrep and 



eqs.(|30|) are the "fourth power" of eqs.(|2£^ 

From (pfl) we deduce that whatever is the number of 8-dimensional irreps, the symmetry of I'Saf) 
depends only by the presence or absence of the 4-dimensional one. More exactly if N/2 is even |$af) 
belongs to the one-dimensional irrep Ai and if N/2 is odd to the one-dimensional irrep Bi. 

As noted at the beginning of this section, the spin projection on the singlet subspace does not alter 
the above quantum numbers. Therefore we conclude that the W — singlet state |$^]^°) has total 
momentum Ktot = (O7 0) and transforms as a d wave of x^ — y^ symmetry if N/2 is odd and as an s 
wave of x'^ + y^ symmetry if N /2 is even. The same quantum numbers were obtained in Ref. [0| for 
the ground state of the Hubbard model in the opposite, strong coupling regime. This coincidence is a 
further consequence of Lieb's theorem. Since the ground state at half filling must be unique, no level 
crossing is allowed for finite U , and the symmetry of the ground state is the same at weak and strong 
coupling. 
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APPENDIX A: CONTRIBUTIONS TO THE W MATRIX FROM FILLED SHELLS 

The A^^— body determinantal wave functions with 6*2 = that one can build using the orbitals with 

/ 2N - 2 \ ^ 
e(fc) < and half of those with e(fc) — are a set of I ,^ , ) elements. Each represents one of the 

degenerate unperturbed {U — 0) ground state configurations at half filling. First-order perturbation 
theory requires the diagonalization of the W matrix over such a basis. 

The diagonal elements of the W matrix are just expectation values over determinants \ka] kfj [ . . .). 
Such an expectation value is a sum over all the possible pairs of the bielectronic elements of W like 

W{a(3,a(3) ^^U{ka,\nr\ka){kp\nr\kp) 

r 

~ ^ N^ N^ ~ N^ ~ N^' ^ ' 

r 

the result is independent of fc^ and fc^. Since in any determinant of the set N'^ /2 plane wave states 
are occupied for each spin, there are N^ /A pairs, and the diagonal elements are all equal to UN"^ /A. 
Thus, the diagonal elements shift all the eigenvalues by this fixed amount. 

The off-diagonal elements of the W matrix between determinants that differ by three or more spin- 
orbitals vanish because W^ is a two-body operator. The off-diagonal elements between determinants 
that differ by one spin-orbital are sum of contributions like W{a(3,^(3) = X^r U {ka\nr\k~j) {ki3\nr\kfj) 
that vanish because of the orthogonality of the plane-wave orbitals. One is left with the off-diagonal 
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elements between determinants that differ by two spin-orbitals, whicli coincide with the corresponding 
bielectronic elements W{a(3,jS) — J2r U {ka\nr\k^) {kj3\nr\ks) . This is just the matrix of W over the 
truncated Hilbert space H spanned by the states of the holes in the half filled shell, ignoring the filled 
ones. We stress that there are N — 1 holes of each spin in Shf, thus H is much smaller than the full 
Hilbert space of the Hubbard Hamiltonian; however, since the number of holes grows linearly with N, 
the problem is still far from trivial. 



APPENDIX B: EIGENVALUES OF THE NUMBER OPERATOR 

Here we prove Equation (|4|). Expanding the determinant according to its definition in terms of the 
totally antisymmetric tensor e, 

i\---i\S\ (Bl) 

we see that the term of maximum order in A is (— A)' '; it arises from the fundamental permutation 
(ii, Z2, . . . ) = (1, 2, . . . ). The {\S\ — l)-th order term in A is the sum of \S\ identical contributions also 
arising from the fundamental permutation. Therefore, it is 

H(-A)l^l-. (B2) 

It is not difficult to see that all the other orders in A yield nothing. At order zero one finds 

il---i\S\ 

Since the exponential is totally symmetric in the permutation of zi . . . i\s\ while e is totally antisym- 
metric the sum vanishes. 

Now we analize the first order term in A. One of its contributions is obtained by picking —A in the 



first factor of eq.(Bl), i.e. by setting ii = 1. This contribution can be written as 

i V ei- ■ g'('=2 + ... + fe|s|-fe,2-...-fc,|s|)r- ,g.x 

■|5|-i) Z^ tij2...,|5|e . ^C4j 



Ar2(|5|-i) 

«2---I|S 



Again, the exponential is symmetric «2 • ■ • i\s\ 8'iid the sum in eq.(B4) vanishes. Clearly, this argument 
applies to all the contributions to first order and to all orders less than \S\ — 1. 
So, we are left with the only two nonvanishing terms. 



detK-A| = (-A)l'5|-i(M-A) (B5) 



that is, eq.(|). 



APPENDIX C: SPIN PROJECTION OPERATORS Ps 

Let us consider a determinantal state with 2n spin-orbitals, half of spin up and half of spin down, 
like |I| 

|*)„ = IV'l ■ • ■ V'n)<T|V'«+l • ■ • V'2n)-<T (CI) 

with {tpi\'4>j) — Sij. Due to the anticommuting property of the fermionic operators, {ip) is separately 
antisymmetric in the indices 1, . . . , n and n + 1, . . . , 2n. \'^)a is no eigenstate of S^ and VS* — 0, . . .n, 
Ps\'^)cr 7^ 0. To build Ps, we take advantage of the fact that different S correspond to different irreps 
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of the Permutation Group of the spins. To this end, let us draw the foUowing table, or Young Tableau 



h 12 




in-S 


in-S+1 




in+S 


jl J2 




jn-S 









(C2) 



where the values of the indices ik and jk must be a permutation of 1, . . . , 2n with the constraints 1) 
ik < jk '^k = 1, ... ,n — 5* and 2) ik < ik+i, jk < jk+i- Now we define Cij as the operator that 
exchanges the spins of the \ipi) and \ipj) states; we associate to the above table the operator 



n+S-l n+S-a n-S-1 n-S-b 

n (1+ E ^^..^.+.-.+i) n (1+ E c',,., 

/3=1 



,.)\l{^-c., 



(C3) 



Q = l 



6=1 



This operator antisymmetrizes the indices on the same column and then symmetrizes those on the same 
row. Then Ps is proportional, up to a normalization factor, to the sum of the operators associated to 
all tables (Cl) that comply with the constraints 1) and 2). 

Example: n = 2 =^ |*)^ = \^-^ T V'2 T V's i ^"4 i) = |V'iV''2)Tl'03V''4)i- In this case the three 
projection operators are 



Ps=2 oc 1 2 3 4 



2 3 



3 4 



(C4) 



1 

3 


2 

4 



+ 


1 
2 


3 

4 



Ps=o oc 



If a table has two indices corresponding to states with equal spin in the same column, the action of 
the associated operator on \^)a will yield zero. This means that in the above example we can omit 
the third table of Ps=i and the second one of Ps=o- We have special intererst in Ps=o, which projects 
on the totally antisymmetric irrep. For all n, it always consists of only one table 



P 



s=o 



1 


2 




n 


n+1 


n + 2 




2n 



(C5) 



Using the explicit form of Pg^o one finds that the singlet projection contains l^)^ and |^)-cr in a 
combination, which is symmetric for even n and antisymmetric if n is odd. 
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1 
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1 
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E 


2 
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